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Abstract: We dehne a resonant collision of order /c(> 1) in a family of four-dimensional 
reversible maps. For any specihed k, the bifurcation scenario is the collection of the dif¬ 
ferent possible types of bifurcation of a symmetric hxed point that may be encountered 
through various choices of parameters describing the family of maps under consideration. 
We adopt a perturbative approach, coupled with numerical iterations around orbits ob¬ 
tained perturbatively, to explore phase space structures in the immediate vicinity of 
a resonant collision and thereby to obtain a description of the possible scenarios for 
different values of k. The phase space structures typically involve bifurcating periodic 
orbits, families of invariant curves, and tori, and present interesting possibiliries, espe¬ 
cially around the ’secondary bifurcations’ of the periodic orbits (see below). Based on 
the results of the perturbative and numerical approach we conjecture that three distinct 
scenarios are involved for the cases /c = 2, 3, 4 respectively, while there exists a fourth 
distinctive scenario common to all /c > 4, and we present what we believe to be a rea¬ 
sonably exhaustive description of these scenarios. The case k = 1 involves bifurcations 
of hxed points rather than of periodic orbits, and has been investigated numerically in 
a previous paper. 
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I. INTRODUCTION. 


Resonances in 4D symplectic and reversible maps occnrring throngh the collision 
of a qnartet of mnltipliers on the nnit circle and their snbseqnent departnre from the 
circle are, in a sense, connterparts of corresponding resonances in 4D Hamiltonian and 
reversible flows. 

In the context of Hamiltonian flows a resonance of this type involves the collision 
on the imaginary axis of a qnartet of eigenvalnes associated with a hxed point, and their 
snbseqnent departnre therefrom, and is termed the 1 : — 1 resonance, where the minns 
sign indicates that the pairs of eigenvalnes involved are of opposite Krein signatnre 
[1], An alternate nomenclatnre is 1:1 non-semisimple resonance (as opposed to the 1:1 
semisimple resonance where the eigenvalnes do not leave the imaginary axis). There 
is a long and distingnished history of investigations on the bifnrcation associated with 
this resonance, originating in the restricted three-body problem ([2], see [3] for details 
and for historical notes), and this bifnrcation has come to be known as the Hamiltonian 
Hopf bifnrcation [3]. It is enconntered in the vicinity of the Ronth’s critical mass ratio 
in the restricted three body problem [4], as well as in a spinning ‘orthogonal’ donble 
pendnlnm [5], and in an atom placed in a rotating electric held [6]. 

Analogons resonance and bifnrcation phenomena are enconntered in reversible 
hows, occnrring in the theory of non-linear oscillations ([7,8] ; see references therein). 
While the Krein theory is not applicable to non-Hamiltonian reversible hows, still these 
hows share many characteristic featnres of Hamiltonian ones close to symmetric hxed 
points and periodic orbits (see [7,8] for details). 

In this context, it is nsefnl to look into analogons resonances in 4D symplectic and 
reversible maps involving a qnartet of mnltipliers of a hxed point on the nnit circle 
where, generically speaking, resnlts resembling those in hows are obtained. However, 
4D maps are, trnly speaking, representatives of dynamics of 6D hows, and this shows 


2 



up in an essential complication in collisions in 4D maps as compared to 4D flows : in 
case of maps one has to distinguish between the rational and irrational collisions. 

While a collision itself is a resonance phenomenon, a rational collision indicates 
the presence of an additional resonance, and may be termed a ‘resonant’ collision to 
distinguish it from an irrational or ‘non-resonant’ collision, the latter being generic in 
the measure-theoretic sense. 

Irrational collisions in 4D reversible maps have been studied in [9,10,11]. In the 
Hamiltonian setting they have been investigated in [12] with rigorous results. Rational 
collisions at on the unit circle, on the other hand, have been analysed in [13] for 
sympletic maps, while a rational collision at for a certain family of reversible maps 
analysed in [14] revealed features quite similar to the symplectic collision (see also 
[15]). Interesting numerical results on bifurcations near rational and irrational collisions 
are to be found in [16]. 

A collision of a quartet of multipliers at (for a hxed point of a symplectic map 
or a symmetric hxed point of a reversible map), with (/> = 27r//c {k positive integer) will 
be termed a resonant collision of order k. 

While a collision of multipliers at an irrational angle is generally accompanied by the 
bifurcation of families of invariant curves, resonant (or rational) collisions are typically 
characterised by bifurcations of periodic orbits. 

In this paper we present perturbation-theoretic and numerical results for resonant 
collisions of order k [k = 2, 3, ■ • ■, 6), from which we shall see that each of the cases 
k = 2,3,4 involves a characteristic bifurcation ‘scenario’ (see below) while the cases 
/c > 4 present a common but distinct scenario. In a sense, the situation resembles 
the bifurcation of periodic orbits in 2D symplectic maps where resonances of order k 
are seen to lead to different scenarios for /c=2,3 and 4 while there occurs a fourth and 
distinct scenario for any /c > 4 [4, 18]. 
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As mentioned above, the bifurcation of period-4 orbits in 4D symplectic and re¬ 
versible maps has been studied in the literature in some detail [13,14], While the 
dynamics in phase space around the periodic orbits was not explored in these works, 
preliminary results in this regard were presented in [15] where families of invariant 
curves around the bifurcating period-4 orbits were investigated. An interesting feature 
in this context is the so called ‘secondary bifurcation’ [13,15] whereby elliptic invariant 
sets involved in a rational angle bifurcation change stability causing a further change 
in local phase space structure. As will be seen in the following, this feature occurs 
commonly in the vicinity of resonant collisions of other orders as well. 

We plan this paper as follows : 

In section II we specify the family of maps to be studied and explain a few introduc¬ 
tory notions. Sections III, IV and V deal with resonant collisions of order 2, 3 and 4 
respectively, the hrst and the last of these being brief partial summaries, with necessary 
elaborations in a few instances, of results of ref. [15]. Section VI includes a few prelim¬ 
inary results on resonant collisions of orders 5 and 6. Section VII contains a synopsis 
and concluding remarks. 

II. RESONANT COLLISIONS IN TWO-PARAMETER FAMILIES OF 4D 
REVERSIBLE MAPS. 

In this paper we consider situations involving two-parameter families of 4D re¬ 
versible maps Ag e such that : 

(i) Ag^e has a symmetric hxed point at the origin for each q, e in some domain D of the 
parameter space (which is 2D, involving parameters q, e) containing the point q = 0 and 
e = 0; 

(ii) Ao,o has eigenvalues (/c=integer> 2) with each eigenvalue repeated twice (we 

call it a resonant collision of order k) or an eigenvalue —1 with multiplicity 4 (we call 
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it a resonant collision of order 2). 

[Remark : For definitions of ‘reversible map’ and ‘symmetric fixed point’, see ref.s 
[9],[10] ; here A denotes the linearisation of A at the fixed point at the origin, the 
relevant parameter valnes being indicated throngh snbscripts.j 

(iii) for any q, e in the domain D, other than (0,0) the linearisation at the origin 
has distinct eigenvalnes. 

(iv) the Jordan block strnctnre of Ao,o contains either a real Jordan block of order 4 
{k = 2) or two complex Jordan blocks each of order 2 {k> 2). 

It will be seen from the following that for each family of maps of the above type, the 
parameter space is typically divided into several regions by two cnrves as shown in 
fig.l (these parameter regions have been marked for fntnre reference as Pla, P16, P2, 
and P3), where the parametrisation has been chosen in snch a way that the cnrves are 
represented by the eqnations 

e = 0, + e = 0 (la, 6), 

these cnrves having a tangency at (0,0). The local strnctnre of the phase space aronnd 
the fixed point is distinct in each of these regions. When any one of the bonndary 
cnrves is crossed at a point other than (0,0) in the parameter space the change in the 
local phase space corresponds to a codimension-1 bifnrcation. The point (0,0) in the 
parameter space, where two bonndary cnrves meet, is degenerate, and the associated 
bifnrcation is of codimension 2. 

The local phase space strnctnre in the vicinity of the fixed point for any Ag ,= is 
essentially described in terms of stability type of the fixed point and the existence and 
stability type of the bifnrcating periodic points (a periodic point is said to bifnrcate from 
the fixed point if it tends to and merges with the fixed point as some particnlar point in 
the parameter space is approached). Additionally, a description of families of invariant 
cnrves and tori aronnd the periodic point may have to be inclnded for a reasonably 
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complete description of the local phase space structure. Depending on the nonlinear 
terms, each family will be found to involve a certain number of types of changes in 
the local phase space structure as the different regions around the point ( 0 , 0 ) in the 
parameter space are approached from one another. The possible different types depend 
on the order of the resonant collision. Thus we get different scenarios of bifurcations, 
each scenario involving its own types of change in the local phase space structure. As 
already stated, the situations with /c = 2, 3 and 4 are found to correspond to three 
distinct scenarios and there is a fourth scenario which is common to the situations 
involving all resonant collisions of order k > 5. We describe below the different scenarios 
obtained from our analysis. However, the description is not complete since each scenario 
typically involves certain degenerate types which have been left outside the purview of 
our investigations. These degenerate types are described by normal forms involving 
higher degree terms compared to the normal forms of the nondegenerate cases. While 
we have not attempted a normal form analysis of the resonances in the present work we 
conjecture that, modulo the degenerate cases, the scenarios indicated above do exhaust 
the typical behaviours of families of 4D reversible maps satisfying criteria (i)-(iv) above. 

III. SECOND ORDER RESONANT COLLISIONS AND PERIOD-2 OR¬ 
BITS 

In order to investigate the second order resonant collision i.e., a collision near the point 
— 1 on the unit circle, we consider the two parameter family of maps described in 
the form of the following fourth order difference equation : 

Xn+2 + Xn-2 “ 2(q - 2){Xn+l + Xn-l) + - 4g + 6 + e)Xn = (dx\ (2), 

where e, q are small bifurcation parameters and 7 ^ 0 is a control parameter. The case 
/? = 0 is degenerate and requires higher degree nonlinear terms for its description. We 
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consider only the non-degenerate case and normalise (3 to \ (3 \= 1 throngh a snitable 
scale transformation. The inclnsion of qnadratic terms in the right hand side of eq. (2) 
does not lead to a distinct bifnrcation scenario and so, for the sake of simplicity, we 
omit snch terms. 

A typical period-2 orbit of the map (2) is of the form 

Xn = {-iTho (3a), 


with 


2 ^ + T 


^>0 = 




(36), 


It follows from eqn.(3b) that a period-2 orbit can bifnrcate from the hxed point in two 
ways, and the corresponding bifnrcation diagrams are presented in hg.s 2-3, each type 
being associated with specihc signatnres of the control parameter (3. 

In each hgnre the horizontal line depicts the variation of q for hxed e, while the 
distance of each point of the bifnrcating branch from the g-axis depicts the amplitnde of 
the orbit for the given q and e. The horizontal q-axis is drawn as a solid (broken) line to 
indicate that the symmetric hxed point abont which the periodic orbit is bifnrcating, is 
linearly stable (nnstable). Similarly, the branch, or portions thereof, denoted with solid 
lines (broken lines) depict linearly stable (nnstable) orbits - anticipating the stability 
resnlts presented below. 

The linear stability analysis of the bifnrcating period-2 orbit is based on the vari¬ 
ational eqnation, obtained by linearising eqn.(2), aronnd a period-2 orbit Xn ■ 

^n-|-2+^n-2 —2((3' —2)(^n-|-l+^n-l) + ('?^ —4(? + 6-|-e)^n = ‘3>(3X‘^^n (^n = X^—Xn) (4). 


In general, for an arbitrary map and for any arbitrary period-two orbit, the variational 
eqnation wonld be a linear diherence eqnation with period-two coefficients, and its 
general solntion wonld be made np of basic solntions of the Floqnet form 


= A"C0 


(8), 
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where A is a quasi-multiplier (see, for an explanation, ref. [17]). In the present case, 
eqn.s 3(a,b),(4) imply that the period-two coefficients are just constants, but we still 
assume a solution of the form (5) to illustrate the procedure in general. Substituting 
eqn.(5) in (4) and equating the coefficients of on both sides, we get 

(A H—-) = —2 q -s/2e (6). 

A 

It is to be noted that eqn.(6) provide us with only one quartet of quasi-multipliers, 
while, in general, there should occur q number of quartets for a q-periodic point of a 4D 
reversible map (see ref. [17]). This is due to the special form of the variational equation 
indicated above, and the second quartet of quasi-multipliers in the present case can be 
obtained just by multiplying the multipliers belonging to the hrst quartet with 

The Floquet multipliers of the period-2 orbits are = A^ {i — 1, ..4) and can, in 
principle, correspond to any one of the four dispositions in the complex plane shown in 
£g.4 (these dispositions are marked R1 through R4 for easy reference). 

Based on these results, we hnd that there arise two distinct bifurcation types {vide 
ref. [15]) which we describe below. 

Type I bifurcation : /I > 0. 

A period-2 orbit exists on the subthreshold side, e < 0, for g > (region Pla 

of parameter space, refer to hg. 1) and q < (-Pl^)? ^md is unstable (disposition 

R3 of hg.4). It merges with the fixed point at e = 0, q = 0 and, on the superthreshold 
side e > 0 (region P3 of hg.l), gets detached from the fixed point with the same type 
of instability {R3). 

Type II bifurcation : /I < 0. 

For e < 0, a period-2 orbit exists if < q < (region P2 of hg.l). From 

eqn.(6) it follows that the period-2 orbit is unstable(with disposition R3 of multipliers) 
for q < — i/—unstable (with disposition R2) for — ^ < q < i/—and stable 



(disposition Rl) for q > y—Further, in this type of bifurcation, there exist no 
period-2 orbits for e > 0. 

The secondary bifurcation 

We see from above that the period-2 orbit appearing in type II bifurcation under¬ 
goes a transition from stable (-R1) to unstable (i?2) conhguration through a secondary 
bifurcation at 



The transition is characterised by the fact that the quasi-multipliers of the period-2 
orbit collide pairwise at ed^=*=®o) when q = go, where 

4sin^^ = ( 8 ). 

This ‘secondary collision’ involves the bifurcation of invariant curves in the vicinity of 
the period-2 orbit. 

Close to the secondary collision, the nonlinear terms become relevant, and the 
deviation = Xn — Xn of a typical orbit from the period-2 orbit under consideration 
is given by the nonlinear difference equation 

^n-r2+^n-2 —2(g —2)(^n-ri+^n-i) + (Q^ —4g-l-6-Fe)^n = + (9). 

and, following [15], we seek solutions of the form 

in = Ar, + + • • ■ ( 10 ). 

Note that the coefficients in eqn.(lO) being period-2 in nature, the coefficients An, Bn, 
Cn, - ■ ■ will also be of period 2, say, 

An = a + A[—1)'^ (llti), 
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B„ =6 + B(-l)" 


(lit), 


C„ = c + C(- 1 )’* 


(11c), 


where a, • ■ ■, (7, • • ■ are constants to be determined. In eqn.(lO), (p is a, rotation angle 
close to the collision angle tt + 9q (this is one of the conjngate angles at which the roots 
Xi (i = 1, ■ ■ ■, 4) collide), say, 

(p^ 71 + 00 + ^ ( 12 ), 

with I <p \« 00. Different valnes of <p wonld lead to a family of invariant cnrves 
organised in islands aronnd the period-2 orbit. 

Snbstitnting eqn.(lO) in eqn.(ll) and eqnating the coefficients of 1, 
terms respectively one gets a system of eqnations in Bn, Cn,--- ■ Snbstitnting 
eqns.(ll a,b,c) in these eqnations and then nsing an order-by-order pertnrbation in 
I q — qo I one hnds that the leading term in Bn is the one involving b and, of all 
the terms in eqn.(lO), this constitntes the dominant contribntion. The next-to-leading 
contribntions in eqn.(lO) are fonnd to arise from terms involving A in An and C in Cn- 
Trnncating the pertnrbation calculation at this order, we get 


A = 


6/3bo 


(q 2 - 2e - 3 q2 


C = 


3/36o 


((q-2 + 2cos2(py - 2e - 3q^) 

^2 {{q ~ 2 — 2cos(pY — 2e — 3^^) 

7 


where 


7 = 


36/3^6o 


21,2 


+ 




(q 2 - 2e - 3 q2) {{q-2 + 2cos2(pY - 2e - 3^2) 


+ 3/3 


(13a), 

(136), 

(13c), 

(14). 


10 




It is found that 7 > 0 and hence there exist two 1-parameter families of invariant curves 
(with associated families of 2-tori) around the period-2 points with rotation numbers 
^ where the two families correspond to (^ > (^1 and ^ < ^2 respectively, with 




2sin9o 


(15). 


The two families merge into a single family as q approaches qq from above and then 


detach from the period-2 orbit as a single one-parameter family of invariant curves 
(together with associated 2-tori) for q < go. In accordance with the results of ref.s 
[9] and [10], the secondary bifurcation can be described as a ‘normal reversible Hopf 
bifurcation’ (an alternative nomenclature is the Naimark-Sacker bifurcation, see [17]). 

Fig.5 shows some of the invariant curves in 2-D projection, obtained on numerical 
iterations, on the (a) subthreshold and (b) superthreshold side (for details, see legend) 
of the secondary bifurcation. These iterations corroborate the results stated above 


relating to the secondary bifurcation associated with second order resonant collisions 
in families of 4D reversible maps. While the invariant curves in hg.(5a) belong to a 
so-called ‘attached’ family, those in £g.(5b) are members of what may be termed a 
‘detached’ family of invariant curves (see sec. 5 below for an explanation). 


IV. THIRD ORDER RESONANT COLLISIONS. 


The 2-parameter family of 4D maps, written in the form of the following 4th order 
difference equation, describes a third order resonant collision at q = e = 0 : 

Xn+2 + Xn-2 “ 2(q - l){xn+i + Xn-i) + {q‘^ - 2q + 3 + e)xn = cxx"^ -|- (3xl_ (16), 

where e and q are small bifurcation parameters and, once again, the origin is a symmetric 
hxed point. 

It is to be noted that the second degree non-linear term in eqn.(16) happens to be 
adequate in describing the bifurcation scenario for the third order resonant collision. 
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since third and higher degree terms are fonnd to lead to no new bifnrcation types while, 
for the second order collision, for which the linearisation of the map at g = e = 0 has a 
distinct Jordan block strnctnre, the third degree term is necessary to describe the fnll 
scenario. Hence we set /? = 0 in the following, thongh in onr nnmerical iterations a 
non-zero valne of jS has been adopted. 

The mnltipliers of the origin nndergo a collision at ^ on the nnit circle when 
q = 0 and e = 0 and, for small q, e, bifnrcating period-3 orbits exist close to the hxed 
point. A typical bifnrcating period-3 orbit is of the form 

Xn — a-\-be 3 +be 3 (17), 

where a and b can be obtained by snbstitning eqn.(17) in eqn.(16) and eqnating the 
coefficients of the constant term and of respectively. 

One can evalnate these constants pertnrbatively in terms of the bifnrcation param¬ 
eters. For the pnrpose of the present paper where we are interested only in describing 
the different bifnrcation scenarios, we retain only the terms of the leading order in the 
bifnrcation parameters. 

Snch a leading order calcnlation shows that, for any given q and e there occnr, apart 
from the trivial solntion a = 6 = 0 , two distinct non-trivial solntions for b, associated 
with a single solntion for a. However, both the solntions are fonnd to correspond to 
the same period-3 orbit when snbstitnted in eqn.(17). A fnrther observation is that the 
bifnrcating orbit is independent of the signatnre of a, and we can chose 

a > 0 

with no loss of generality. In snmmary, there occnrs only one bifnrcating period-3 orbit 
for each point in parameter space, and the valnes of a, 6 , describing this orbit are (with 
a chosen positive) 
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b = b 


(18a), 



A = (—4 + 6q^ — 6e) + {6eq — 2q^) + -(e + (20a), 

3 

291 

B = {4- 12q^ + 12e) + (4^^ - 12eQ) - -^(e + q'^f (206). 

9 

Thus, the multipliers of the period-3 orbit are obtained from the expression 


A-I- — = 2-l-3e — 3q^ ± bA/e^ + q^ + e q'^ 


It is easily seen from eqn.(21) that, for arbitrary choice of the parameters e and one 
pair of multipliers of the period-3 orbit lies on the unit circle and the other on the real 
line, corresponding to the disposition R3 in hg.4. In other words, the bifurcating orbit 
is unstable in its domain of existence. 

The bifurcation diagram following from these results is presented in hg.6. One hnds 
that, for e < 0 (the subthreshold side of the bifurcation), the period-3 orbit merges 
with the hxed point for | q'^ + e |—^ 0 while, for e > 0 (superthreshold side) the orbit 
is bounded away from the hxed point, the orbit being unstable in both the cases. This 
diagram corresponds to what has been termed a bifurcation of type III in the context 
of fourth order resonant collisions in ref.[15] (see next section). 
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In summary, while the bifurcation scenario for the second order resonant collision 
consists of two types of bifurcation (type I and type II of sec. Ill), only one type (analo¬ 
gous to type III of ref. [15]) is found to occur in the bifurcation scenario associated with 
the third order collision. In hg.7 we show the 2D projection of typical trajectories start¬ 
ing close to the period-3 orbit near the bifurcation, where one hnds that the trajectories 
do diverge away from the orbit, in conformity with its unstable nature. 

V. RESONANT COLLISIONS OF ORDER FOUR 

In order to describe the bifurcation of period-4 orbits arising from a collision of 
multipliers at (/c = 4), we consider, as in ref.[15], the following two-parameter family 
of maps Ag e, written in the form of a 4th order difference equation : 

Xn+2 + Xn-2 + MiXn+l + X^-l) + {4q‘^ + 2 + = P(xl -fXn+lXnXn-l) ( 22 ), 

where q, e are small bifurcation parameters and {3 and 7 are control parameters which 
will be found in the sequel to govern the nature or the type of the bifurcation. These 
parameters will be assumed to satisfy the nondegeneracy conditions 

(3^0 and 7 ^ 7 ^ 1 (23). 

Quadratic terms in eqn.( 22 ) have not been included since they are found not to affect 
the scenario of bifurcations involved. The origin is once again a symmetric hxed point 
of and the multipliers at the origin undergo a collision at at g = 0 and e= 0 . 

A typical bifurcating period-4 orbit of ( 22 ) is of the form 

Xn = (iTbo + i-iTb* (24), 

where bo is obtained by substituting eqn (24) in (22) and equating the coefficients of 
(i)^ and Two distinct types of bifurcating period-4 orbits are obtained in this 

way, which we denote by 0\ and 02 respectively. 

01 : bo = bl = pi (25a), 
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with 


(256). 


2 _ e + 

/? 

02 : bo = p 2 {l + i) (26a), 

2 

with pI = (266). 

As explained in [15], the two types of orbits have distinct symmetry properties nnder 
the reversing involntions pertaining to ( 22 ). 

The linear stability analysis for the bifnrcating period-4 orbits can be done in 
the same way as that for period -2 orbit from the appropriate variational eqnations. 
This has been presented in detail in ref. [15] and will not be repeated here. When 
the stability resnlts of [15] are combined with eqn.s (25a) throngh (26b), one hnds 
that the bifnrcation scenario associated with the fonrth order resonance involves three 
bifnrcation types. We designate these as bifnrcations of type I, II, and III, corresponding 
respectively to bifnrcation diagrams of hg.s 8 , 9, and 10. A brief description of these 
types, corresponding to different sets of valnes of the control parameters (3 and 7 , is as 
follows. 

Type I : /I > 0, 7 < 1 

For 7 < —1, on the snbthreshold side (e < 0), both the bifnrcating branches of 
period-4 orbits (01 and 02) exist for q > and q < (i.e., in regions Pla 

and P16 of the parameter space, refer to hg.l). In this sitnation 01 is linearly stable 
(disposition R1 of hg. 4) and 02 is nnstable (disposition i?3) (hg.Sa). For e = 0, they 
merge together at the origin at g = 0 (hg. 8 b) and for e > 0 , they get detached from the 
origin with their respective stability types remaining nnchanged (hg. 8 c). The sitnation 
remains essentially the same for — 1 < 7 < 1 , with only the stability types getting 
interchanged (disposition R?> for 01 and R1 for 02). 

Type II : /I < 0, 7 < 1 

For e < 0, two branches exist simnltaneonsly for < q < and are 
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globally connected (as the second bifnrcation parameter q is varied), getting annihilated 
simnltaneonsly as e becomes positive. 

Here we enconnter once again the interesting sitnation that one of the two branches 
nndergoes a change of stability at some critical valne of q (for hxed e < 0 ) throngh 
a collision of Floqnet mnltipliers on the nnit circle near +1 (Fig.9). Thns, for 7 < 0, 
the branch 02 is nnstable (disposition R3) thronghont the range —\/\ e | < q < a/| e | 
while, for the 0\ solntion, there occnrs a transition from linearly stable (disposition 
i?l) to nnstable(disposition R2) conhgnration as q‘^ crosses the threshold 


1 - 7 + 2a/2v/-1 - 7 
^5-7-2v^v/=r^ 


(27) 


from above, dne to a ‘small angle’ collision (see [19] for an explanation) of Floqnet 
mnltipliers. As q‘^ is made to decrease farther, a second transition occnrs at 


2 2 1 - 7-2^2^-! -7 

Q = Q2 = - r- / : 

5-7-2y2V-l - 7 


(28), 


the disposition of the mntipliers changing over to R4 which, however, does not change 
the nnstable character of the branch (this transition has not been shown in £g. 9). For 
0 < 7 < 1, on the other hand, the 01 orbit is nnstable (disposition R3) thronghont 
the range —a/| e | < q < a/| e |, while for the 02 branch there occnrs a transition from 
linearly stable(i?l) to nnstable {R2) conhgnration as q^ crosses the valne 

1 + 


Qs = -e 


2 - 7 + a/1 - 7 ^ 


(29), 


from above, dne to a small angle collision of the mnltipliers. As q^ is made to decrease 
farther, a second transition occnrs as in the case of the 0\ branch above. 

Type III : 7 > 1. 

In this case one branch {01 for (3 > 0) exists for q > and q < 

(regions Pla and Plb of hg. 1) and the other branch (02) exists for < q < 
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(Fig.10a). Each of the branches is unstable (disposition R3) in its domain of existence. 
For e > 0 the branch 0\ gets detached from the origin and 02 is annihilated (Fig.10c). 
The branches get interchanged if /I < 0. 

Thus, only for the type II collision near (i.e., for /I < 0 and 7 < 1 ), there occurs 
a transition of a period-4 branch from linearly stable (disposition Rl) to an un stable 
(i?2) conhguration through a small angle secondary collision. This is the co-dimension 
1 collision discussed in ref.s [9,10]) and so the full co-dimension 2 collision studied in 
[19] is not relevant here. The transition in the disposition of multipliers from R2 to R4 
referred to above again involves a co-dimension 1 collision, and will not be discussed 
here (see, e.g., [18], section 36). 

The phase space dynamics in the vicinity of the secondary collision can once again 
be investigated by going beyond the linear terms. The analysis is essentially similar to 
the one outlined in case of the second order resonant collision and has been presented in 
detail in ref. [15]. The principal result in this respect is that the secondary bifurcation 
is superthreshold or ‘normal’ in character (see [9,10] for an explanation). 

One can proceed further from this result and construct perturbatively the ‘islands’ 
of invariant curves surrounding the period-4 points in the vicinity of the secondary 
collision. Thus, on the subthreshold side of the bifurcation, each linearly stable period- 
4 point is surrounded by a family of invariant curves (forming an island, there being 
four such islands for the period-4 orbit under consideration), members of the family 
passing arbitrarily close to the period-4 point for rotation angles chosen sufficiently 
close to This family persists on the superthreshold side of the bifurcation, but now 
detached from the period-4 point under consideration, i.e., there is now a lower bound 
to the distance of orbits belonging to this family from the period-4 point. These two 
types of families, with zero and non-zero lower bounds to distances from the hxed point 
or periodic point under consideration will be termed ‘attached’ and ‘detached’ families 
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respectively. 

Attached families of invariant cnrves are commonly fonnd to occnr in reversible 
mappings aronnd linearly stable fixed or periodic points (see ref. [15]), while there 
exists no standard resnlt in the literatnre concerning detached families. As described 
in ref.s [10,11,15], snch detached families are to be observed in the vicinities of ‘nor¬ 
mal’ reversible Hopf bifnrcations. We have already come across attached and detached 
families of invariant cnrves aronnd period-2 points in sec. III. 

Fignres 11 and 12 represent resnlts of n um erical iterations with the family of maps 

(eqn. (22)) close to the secondary bifnrcation. Fignre 11a shows three invariant 
cnrves belonging to an attached family forming an island aronnd one of the period-4 
points of 01 on the snbthreshold side of the bifnrcation (for details, see legend) in a two- 
dimensional projection. Fignre 11(b) on the other hand shows three invariant cnrves 
belonging to a detached family on the snperthreshold side of the secondary bifnrcation 
aronnd the 01 orbit. Fignres 12(a,b) show members belonging to corresponding families, 
on the snbthreshold and snperthreshold sides respectively, for the 02 branch (again,see 
legend, for details of each hgnre). 

Two-dimensional projections of trajectories initiated very close to the two types of 
period-4 orbits (01 and 02), on snbthreshold and snperthreshold sides of each of the 
three types of bifnrcation have been presentend in hg.s 13, 14, and 15 respectively. In 
each hgnre the persistence of the trajectory signihes stability of that orbit, whereas a 
precipitate departnre away from the orbit signihes its instability. Thns, hg. 13(a,b) 
represent 2D projections of trajectories initiated very close to the two types of period-4 
orbits in type I bifnrcation, where one hnds that the trajectory aronnd 01 persists near 
it forming an island, while one initiated near 02 departs away from it, following a path 
close to the separatrix. 

Similarly, hg. 14(a-b) depict trajectories initiated near the 01 and 02 orbits in type 
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II bifurcation, now close to the secondary bifurcation. While hg. 14(a) corresponds to 
the subthreshold side of the secondary bifurcation, in which the trajectory initiated close 
to 01 forms an island and the one initiated close to 02 makes a circuit in the vicinity 
of the separatrix, hg. 14(b) corresponds to the superthreshold side, where trajectories 
are seen to depart away from both 01 and 02 orbits. 

Fig. 15 presents similar trajectories for type III bifurcation in the fourth order 
resonant collision. As seen from hg. 10, there exists only one type of orbit {01 or 
02) for each specihcation of parameter values {q, e) on either side of the bifurcation 
(e < 0, e > 0), and that orbit is unstable. Fig. 15(a, b) show trajectories initiated close 
to the 02 and 01 orbits respectively for e < 0, and for two distinct values of q, and in 
each case the trajectory is seen to diverge away from the period-4 orbit concerned. Fig. 
15(c), on the other hand depicts a similar trajectory for e > 0, where only the 01 orbit, 
now detached from the hxed point, survives. 

VI. HIGHER ORDER RESONANT COLLISIONS. 


Collisions of order five. 

We consider hfth order resonant collisions with reference to a family of maps 
written in the form of a fourth order diherence equation. 


Xn+2 + Xn-2 — 4cos())o(-|- Xn-i) -|- {Acos'^fio + 2 + e)Xn — cxX^ -|- fiX^ (30)? 


where now cosfio = cos^ + '?/2, e and q being small bifurcation parameters as in the 
previous sections. 

The multipliers of the origin undergo a collision at at q = 0 and e = 0. A 
typical bifurcating period-5 orbit is of the form 


_ 2n7ri ^ — 2n7ri 4n7vi ^ —4n,7ri 

Xn — a + be ^ -{-be 5 -l-ces -{- c e 5 


(31). 
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ei ^ e + 


€2 = 6 + 4{cos(f)o ~ cos — 

5 


(346), 

(34c), 


respectively. 

These equations tell us that, in the leading order of perturbation, there is only one 
period hve orbit for any given value of the parameters e and q. However, on going over 
to the next order of perturbation (i.e., to terms of the second degree in e and q) one 
hnds that there exist, in fact, two distinct types of period-5 orbits, which we designate 
as 01 and 02 respectively, distinguished by the value of the coefficient 6 in eq. (31) : 
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01 : 

with 
02 : 
with 

Where pu {i = a, b) is given by 


b = b*= Pa 
Pa = PO + Pla 
b = pbe 5 

Pb = Po + Plb 


with 


Pii 


plxcos4:(f)i 

2ecos4>i + Axp^cosA(^i 


[i — a, b) 


(35a), 

(356). 

(зба) , 

( збб ) . 


(37), 


= 0 , (pb 


7T 

5 


38(a, 6), 


and 


X 


3a3 

— H- 

€2 


(39). 


In other words, there is a degeneracy between the two types of period-5 orbits in the 
leading order of pertnrbation, which is removed in the next higher order. It is also 
apparent from eqn. (39) that the qnadratic term in eqn. (30) is important in this 
removal of degeneracy. 

Having obtained the pertnrbative expressions for the period-5 orbits it is a rontine, 
if tedions, affair to follow trajectories initiated close to the orbits for varions sets of 
parameter valnes, scanning across the parameter space. Pertnrbation calcnlations of 
the type resorted to in the case of lower order resonant collisions present increasing 
magnitndes of difhcnlty as one comes across resonant collisions of order 5 and above. 
Hence, in absence of a pertnrbative stability analysis of the orbits and a pertnrbative 
constrnction of the families of invariant cnrves close to these orbits, we invoke the 
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approach of numerical construction of trajectories close to the periodic orbits so as to 
have an idea of the bifurcation scenarios involved. 

In the present instance, such an approach does seem to lead to a reliable picture, 
which we summarise as follows : 

Period-5 orbits can bifurcate from a hxed point in a hfth order resonant collision in 
either of two ways that resemble the type I and type II bifurcations in a fourth order 
resonant collision. Thus, in type I bifurcation there exists a stable and an unstable orbit 
on either side of the bifurcation (e < 0, and e > 0), while in type II bifurcation, a pair 
of orbits exists only on the subthreshold side e < 0, there being no period-5 orbit on 
the superthreshold side. Of the members of the pair bifurcating from the hxed point for 
e < 0, one is always unstable, while the other undergoes a secondary bifurcation which 
is ‘normal’ in nature (see sec. V). 

Fig. 16 presents trajectories initiated close to the period-5 orbits in a type I bi¬ 
furcation for e < 0 (hg. 16a) and e > 0 (hg. 16b). One does hud from the hgure that 
the trajectory near one of the orbits belongs to an island of invariant curves, signifying 
that the orbit concerned is stable, while the trajectory near the other orbit moves away 
following a separatrix loop, and that essentially the same picture persists on the two 
sides of the bifurcation. 

Fig. 17, on the other hand, depicts a secondary bifurcation associated with a 
type II bifurcation where one hnds that trajectories on the subthreshold side of the 
secondary bifurcation (hg. 17a) are similar to those in the type I bifurcation while, 
on the superthreshold side (hg. 17b), trajectories move away from both the orbits, 
signifying that they are now unstable. 

Collisions of order six. 

Bifurcations of period-6 orbits can be studied with the family of maps (^g,e) written 
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in the form 


Xn+2 + Xn-2 — 2(1 + q){XnJ^i + X^-l) + + 2(J + 3 + e)Xn — aX^ + (40), 


where, as before, e and q are small bifnrcation parameters, and a, (3 are control param¬ 
eters. 

As seen from eq. (40) the mnltipliers of the origin nndergo a collision at at 
q = 0 and e = 0. A typical bifnrcating period-6 orbit is of the form 



.2n7ri _ -2n7ri 4mri ... -4mfi 

Xn — a + d{—l) +he e. -\-be s qq a c e s 

( 41 ), 

where a. 

d, b, b*, c, c* can be obtained pertnrbatively. In the leading order of pertnr- 

bation in 

e and q^ one hnds 

2a\ b ^ 

ao ^ 

Co 

(42a), 


ab'^ 

Co ~- 

62 

(426), 


Co ^ 

62 

(42c), 


5(b^ + b*^) 

do ^ 

63 

(42d), 

and 

1 ^1= Po 

(42e), 

with 

2 

Po = - 

(43a), 

where 

y = 2ct^(-1-) -|- 3/3 

60 62 

(436), 


60 = 6 + (<? - 1)^ 

(43c), 


ei ^ e + q'^ 

(43d), 
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(43e), 


e 2 — e-\- {q-\- 2)^ 

€3 = e + (q + 3)^ (43/), 

2ct^ 

5= -+/3 (43<7). 

£2 

Thus, there exists only one bifurcating orbit in the leading order of perturbation 
which, however is to be interpreted as just a degeneracy between distinct orbits that is 
lifted in the higher orders. Indeed, in the next order of perturbation, one hnds that as 
in the case of a resonant collision of order hve there exist two distinct types of period-6 
orbits, which we describe as follows : 


01 : 


b = b*= Pa 


(44a), 


with 
02 : 
with 


Pa — Po + Pla 
b = pbee 

Pb = Po + Plb 


Pli {i — a, b) being given by 


Pli 


Po(xicos/>i -^X2Cos5(/j) 
2ecos(t)i 5po(xicos(/i -|- X2Cos5(j)i 


where 


= 0 , 


= TT 


Xi = 


6 


2a‘^6 3(35 6a‘^(3 12a‘^(3 12a‘^(3 


+ 

£2^3 ^3 


+ 




+ 


+ 


£0^2 


2a‘^d 3(35 3a^(3 

X 2 =-^-^-^ 


(446), 

(45a), 

(456), 


(46), 

(47a, 6), 
(47c), 

(47d). 


£263 £3 ei 

Once again, having obtained the perturbative expressions for the period-6 orbits, the 
dynamics in phase space close to these orbits can be studied through numerical iterations 
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in order to construct the bifurcation scenario. The result emerging from such an exercise 
is : 

The scenario in the sixth order resonant collision is analogous to that in the fifth order 
collision, consisting of two types of bifurcation resembling in turn the type I and type II 
bifurcations associated with the fourth order collision. In particular, the type II bifur¬ 
cation involves a secondary collision whereby one of the two period-6 orbits undergoes a 
transition from the linear stablility to instability. 

Evidence for the two types of bifurcations in the family of maps given by eqn.(40) 
is presented in the form of hg.s 18 and 19, corresponding in signihcance to hg.s 16 and 
17 respectively for the two types of bifurcation in the hfth order resonant collision (see 
legends for details). 

These observations based on numerical iterations, together with the results of the 
previous sections lead to a conjecture concerning resonant collisions in families of re¬ 
versible mappings presented in the next section. 

VII. SUMMARY : A CONJECTURE. 

In this paper we have considered resonant collisions of the multipliers at a sym¬ 
metric hxed point of families of 4D reversible maps specihed as in sec. II. A resonant 
collision has been dehned in sec. I as a collision of multipliers on the unit circle at angles 
positive integer). The case k = 1 involves the bifurcation of fixed points, as 
distinct from the cases k > 1 where periodic orbits are found to bifurcate from the 
hxed point under consideration, and has been investigated from a numerical point of 
view in ref. [19]. The present paper addresses resonant collisions with k > 1. The 

families of maps we have considered for different specihed values of k, while conforming 
to the criteria mentioned in sec. II, have all been chosen to be of the so-called de Vo- 
gelaere type (see, e.g., ref. [9]), these being relatively easy to handle analytically and 
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numerically (eqn.22, corresponding to k = 4, is an exception since the de Vogelaere 
form does not lead to the complete scenario in this case). 

Collecting the results obtained from perturbative calculations and numerical iterations 
we formulate the following conjecture : 

There exists a characteristic bifurcation scenario for each of the cases k = 2, 3,4, and a 
fourth and distinct scenario for k > A. 

The scenario for k = 2 involves two bifurcation types as depicted in fig.s 2 and 3, 
and in each of these there exists only one bifurcating period-2 orbit. This orbit undergoes 
a secondary bifurcation at a particular value of q (0 < q < — e) ), with attached and 

detached families of invariant curves on the subthreshold and superthreshold sides of the 
bifurcation respectively. 

There is only one bifurcation type involved in the scenario for k — 3 (fig. 6), the 
associated period-3 orbit being unstable on either side of the bifurcation. This bifurcation 
is analogous to one of the three types involved in the scenario for k = f. 

Three bifurcation types are involved in the scenario for k — 4 (fid-s 8, 9, 10). We 
have termed these type I, type II, and type III bifurcations respectively. In one of these 
(type II) there occurs a secondary bifurcation of one of the two associated period-f orbits 
at two symmetrically located values of q for any given e < 0. 

The scenario for any k > 4 involves two bifurcation types analogous respectively to 
types I and II for k = 4 (fig.s 8, 9). 

The bifurcating periodic orbits are associated with families of invariant curves which 
may be of either ‘attached’ or ‘detached’ type (see sec. V for explanation) corresponding 
to whether the orbit under consideration is linearly stable or is an unstable orbit on 
the superthreshold side of a secondary bifurcation. Our perturbative and numerical 
approach gives us a good idea of these families sufficiently close to the resonant collisions. 

Additionally, one can talk of 2-dim tori in the immediate neighbourhoods of the 
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bifurcating orbits. Analogous to the families of invariant curves, these tori can also 
be of the ‘attached’ and ‘detached’ types. Thus, each of hg.s 20(a,b) presents a two- 
dimensional projection of a torus around a period-4 orbit, respectively on the subthresh¬ 
old and superthreshold sides of the secondary bifurcation in a resonant collision leading 
to a type II bifurcation. 

In other words, it is possible to have considerable information of the dynamics 
in phase space close to a resonant collision of multipliers of a symmetric hxed point 
of a family of 4D reversible mappings from the perturbative approach coupled with 
numerical iterations of the type presented in this paper. One hopes that this will help 
in the formulation of rigorous results in this largely uncharted area [20]. 

Acknowledgement. We thank M. B. Sevryuk for suggesting the present area of in¬ 
vestigation to one of us (A.L). 
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FIGURE CAPTIONS. 


Fig. 1 Regions in the (e — q) parameter sapce divided by the cnrves (t) e = 0 and 
(ii) q‘^ + e = 0. 

Fig. 2 Schematic bifnrcation diagram for Type I bifnrcation {(3 > 0) of a period-2 
orbit aronnd a symmetric hxed point for (a) e < 0, (b) e = 0, (c) e > 0. For explanation, 
see text, Sec. III. 

Fig. 3 Schematic bifnrcation diagram for Type II bifnrcation {(3 < 0) of a period-2 
orbit aronnd a symmetric hxed point for e < 0. For explanation, see text. Note the 
occnrence of the secondary bifnrcation. The snperthreshold side (e > 0) is not shown 
in the hgnre since there exists no period-2 orbit on this side. 

Fig. 4 Different possible disposition of the Floqnet mnltipliers at a symmetric hxed 
point or periodic orbit for 4D reversible maps. For explanation, see text. Sec III. 

Fig. 5 (a) 2-dim projection of a set of invariant cnrves aronnd one of the period-2 

points in second order resonant collision (type II bifnrcation) on the snbthreshold side of 
the secondary bifnrcation, with e = —10“^, 13 = —1 and q = 2.585 x 10“^, corresponding 
to (^ = 2.17 X 10“^, 2.37 X 10“^ and 2.57 x 10“^ respectively. (b)Projection of a set of 
invariant cnrves aronnd the period-2 point on the snperthreshold side of the secondary 
bifnrcation, with q = —2.58 x 10 ^ for 1 X 10 and 8 X 10 ^ respectively, other 
parameters being the same as in (a). 

Fig. 6 Schematic bifnrcation diagram of period-3 bifnrcation 5 (a) c ^ 0, (b) € — 0, 
(c) e > 0. 

Fig. 7 2-dim projection of a trajectory initiated close to a period-3 point near the 
resonant collision with (3 — a = 1 corresponding to (a) e = —1 x 10“^, q = 8 x 10“^ 
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and 5 = 1 X 10“^, and (b) e = 1 x 10“®, q — 8 x 10“^ and 5 = 1 x 10“^; here and in 
the snbseqnent hgnres the parameter 6 specihes the initial condition for the trajectory; 
we hrst calculate period-3 orbit (Xn) from eqn.(17) and then take the initial condition 
of the trajectory as = Xn -1- n(—(n = 1, 2, 3, 4). 

Fig. 8 Schematic bifurcation diagram for Type I bifurcation (/9 > 0, 7 < 1) of 
period-4 orbits around a symmetric hxed point for (a) e < 0, (b) e = 0, (c) e > 0. The 
labels 0\ and 02 are shown for 7 < 0 while for 0 < 7 < 1 they should be interchanged. 
For explanation, see text, sec.VI. 

Fig. 9 Schematic bifurcation diagram for Type II bifurcation {(3 < 0, 7 < 1). Only 
the subthreshold situation (e < 0) is shown. The orbits shrink to the hxed point at 
e = 0, and cease to exist for e > 0. The labels 01 and 02 are shown for 7 < 0 while for 
0 < 7 < 1 they should be interchanged. 

Fig. 10 Schematic bifurcation diagram of Type-Ill bifurcation (7 > 1); (a) e < 0, (b) 
e = 0, (c) e > 0. The labels 01 and 02 are shown for (3 > 0 while for P < 0 they should 
be interchanged; see text, sec.VI. 

Fig. 11 (a) 2-dim projection of a set of invariant curves around one of the period-4 

points of the 01 orbit on the subthreshold side of the secondary bifurcation, with e = 
— 10“^, P — —1, 7 = —3.02 and q = —2.585 x 10“^ corresponding to p — 1.173 x 10“^, 
1.178 X 10-4 and 1.183 X 10 4 respectively, (b) Projection of a set of invariant curves 
around one of the period-4 points of the 01 orbit on the superthreshold side of the 
secondary bifurcation, with q = —2.5835 x 10“^ for = 1 x 10“®, 3 x 10“^ and 5 x 10“^ 
respectively, other parameters being the same as in (a). 

Fig. 12 (a) 2 -dim projection of a set of invariant curves around one of the period- 

4 points of the 02 orbit on the subthreshold side of the secondary bifurcation, with 
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e = —10 /? = —1, 7 = —10 ^ and g =—2.578x10 ^ corresponding to (^ = 6.99x 10 

7.02 X 10“® and 7.05 X 10 ^ respectively, (b) Projection of a set of invariant cnrves 
aronnd one of the period-4 points of the 02 orbit on the snperthreshold side of the 
secondary bifnrcation, with q = —2.577 x 10“^ for (^ = 1 x 10“®, 2.2 x 10“® and 
3.4 X 10“^ respectively, other parameters being the same as in (a). 

Fig. 13 2-dim projection of two trajectories starting from close to the period-4 orbits 
(01 and 02) in type I bifnrcation near the resonant collision with /? = 1, 7 = 0 and 
q = 2 X 10“^, corresponding to (a) e = —1 x 10“^, h = 1 x 10“^ for 02 and h = 0 for 
01 orbit and (b) e = 1 x 10“^, h = 1 x 10“^ for 02 and h = 0 for 01 orbit; we hrst 
calcnlate period-4 orbit {X^) from eqn.(24) and then take the initial condition of the 
trajectory as = X^ -|- h (h being taken nonzero only for n = 4). 

Fig. 14 2-dim projection of two trajectories starting from close to the period-4 orbits 
(01 and 02) in type II bifnrcation near the secondary bifnrcation with /? = —1, 7 = 0 
and e = —1 X 10“^, corresponding to (a) q = 9.5 x 10“^, h = 5 x 10“® for 02 and 
h = 1 X 10“® for 01 orbit and (b) q = 8.1 x 10“^, 5 = 1 x 10“® for 02 and h = 1 x 10“® 
for 01 orbit ; we hrst calculate period-4 orbit {X^) from eqn.(24) and then take the 
initial condition of the trajectory as X^ = X^ + S {6 being taken nonzero only for n = 4). 

Fig. 15 2-dim projection of a trajectory starting from close to the period-4 orbit (01 
or 02) in type III bifurcation near the bifurcation point with /9 = 1, 7 = 2 corresponding 
to (a) e = -1 X 10-4, q = -9x IQ-^ and h = 2.5 x IQ-^ for 02, (b) e = -1 x IQ-^, 
q — —1.2 X 10“^ and 5 = —1 x 10“® for 01 and (c) e = 1 x 10“^, q = —1.2 x 10“^ and 
h = —1 X 10“® for 01 orbit; we hrst calculate period-4 orbit (Xn) from eqn.(24) and 
then take the initial condition of the trajectory as Xn = Xn + n{—l)'^dcos{nj). 

Fig. 16 2-dim projection of two trajectories starting from close to the period-5 orbits 
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(01 and 02) in type I bifnrcation near the bifnrcation point with /3 = —.5, a = 1 and 
q = 6 X 10“^ corresponding to (a) e = —1 x 10“^ h = 1 x 10“^ for 02 and h = 0 for 
01 orbit and (b) e = 1 x 10“^, h = 5 x 10“^ for 02 and h = 1 x 10“® for 01 orbit; we 
hrst calcnlate period-5 orbit (Xn) from eqn.(31) and then take the initial condition of 
the trajectory as = X^ -|- h (h being taken nonzero only for n = 4). 

Fig.17 2-dim projection of two trajectories starting from close to the period-5 orbits 
{Ol and 02) in type II bifnrcation near the secondary bifnrcation with jS = —2, a = 1 
and e = —5 X10“^ corresponding to (a) g = 6 x 10“^ h = 5 x 10“^ for 02 and h = 1 x 10“^ 
for 0\ orbit and (b) q = 4.806 x 10“^, h = 0 for 02 and 0\ orbit ; we hrst calcnlate 
period-5 orbit (X^) from eqn.(31) and then take the initial condition of the trajectory 
as Xn = Xn -|- h (h being taken nonzero only for n = 4). 

Fig. 18 2 -dim projection of two trajectories starting from close to the period -6 orbits 
(Ol and 02) in type I bifnrcation near the bifnrcation point with (3 = 1, a = 1 and 
q = —2 X 10“^ corresponding to (a) e = —1 x 10“^ h = 1 x 10“^ for Ol and h = 0 for 
02 orbit and (b) e = 1 x 10“^, 6 = 2 x 10“^ for Ol and h = 0 for 02 orbit; we hrst 
calcnlate period-5 orbit (X^) from eqn.(41) and then take the initial condition of the 
trajectory as Xn = Xn + 2n{—l)'^dcos{2n4>) {4>=0 for Ol and ^ for 02). 

Fig. 19 2 -dim projection of two trajectories starting from close to the period -6 orbits 
(Ol and 02) in type II bifnrcation near the secondary bifnrcation with jS — —4, a = 1 
and e = —1 X 10“^ corresponding to (a) q = — 6.8 x 10“^ S — —1.6 x 10“® for 02 and 
h = 1 X 10“® for Ol orbit and (b) q — —6.2 x 10“^, h = 0 for 02 and Ol orbit ; we 
hrst calcnlate period -6 orbit (Xn) from eqn.(41) and then take the initial condition of 
the trajectory as Xn = Xn + 2n{—l)^dcos{2n4>). {4>=0 for Ol and ^ for 02). The 
inset shows one of the stable period -6 orbit snrronnded by the sepratrix of the un stable 
period -6 point. 
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Fig. 20 (a) 2-dim projection of a trajectory around one of the period-4 points of 

the 0 \ orbit on the subthreshold side of the secondary bifurcation, with e = — 10 “^, 
(3 — —1, 7 = —3.02 and q = —2.585 x 10“^ ; we hrst calculate an initial condition 
(Xi, X 2 , Xs, X 4 ) corresponding to one member belonging to a family of invariant curves, 
with (j) = 1.18 X 10“^ (see ref.s [11,15] for explanation ); a slightly different initial 
condition is then taken with Xi, X 2 , X^ unchanged and X 4 replaced by X 4 - 1 - h (h = 
3 X 10“^®) and iterations performed; the trajectory winds on a 2-torus around the 
invariant curve, (b) A 2-torus on the superthreshold side near a period-4 point of the 
01 orbit ; with e = —10“^, [3 — —1, 7 = —10“^ and q = —2.577 x 10“^, ^ = 10“^, 
d = 10 “^ (see legend for (a) ). 
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